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ABSTRACT 


The Electrocardiogram (ECG) record is a procedural electrical activity of the heart which is noninvasive recording and is acquired by 
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surface electrodes at designated locations on the skin of patient's body during which noise such as power-line interference (PLI) with 
60 Hz frequency is acquired from power mains. In this research paper a new ECG enhancement design using multiband structured 
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sub band adaptive filter (MSAF) is constructed to solve structured problems in conventional sub band adaptive filters. This paper 
investigates the new detailed ANC system for ECG signals with robustness based on UFB and NUFB structured MSAF using LMS 
algorithm. The proposed model is potentially a new realization form of algorithm which guaranteed a more stable transformation. 


The goal of the research is to provide solution in order to enhance the performance of ANC in terms of filter parameters which are 
acquired with the help of UFB &NUFB structured MSAF’s using LMS algorithm. Computer simulation demonstrates that the 
proposed system gives improved performance and achieves good adaptation. NUFB structured MSAF algorithms are applied on ECG 
records obtained from standard MIT-BIH data base and the performance is compared with UFB structured MSAF algorithms in terms 
of parameters SNR, MSE, RMSE and distortion. The SNR for various NUFB structured MSAF’s was found to be higher than the UFB 
structured MSAF’s. The five channel NUFB structures with decimation factors (16,16,8,4,2) has on average SNR of 21.7097dB is 
obtained using LMS algorithm which is superior to existing algorithms. 


Keywords: LMS, FB, UFB, NUFB, SAF, MSAF 


1. BACKGROUND 


Signals play essential part in the subject of electrical, digital, communication engineering and medical. The signals associated with 
scientific field are known as biomedical signals. Distinct biomedical alerts like Electrocardiogram, Encephalogram and 
Electromyogram are used for diagnosis as they contain lot of facts. ECG is a critical biomedical instrument for the diagnosis of 
cardiovascular problem that suggested to the electric interest of coronary heart recorded by way of skin electrode. The ECG is a 
bioelectric signal, which statistics the electrical pastime of coronary heart as opposed to time. Therefore, it is a vital diagnostic device 
for assessing coronary heart function [1,2]. The ECG record represents a very necessary live utilized by doctors, because it provides 
very important data a few patients’ viscous condition and general health. The graphical record signal delivers data of the human 
heart, heart position, impulse starting and propagation, heart pattern and transmission of interference and drug effect on 
heart. Usually the band of the graphical record signal is zero.05 rate toa hundred rate. The goal of any ffilter is to 
extract helpful data from strident information. A normal fixed filter is designed in advance with knowledge of the statistics of both 
signal and the unwanted noise but if the statistics of the noise are not known priori, or amendment over time, the coefficients of the 
filter cannot be specified in advance. In these things, adaptive algorithms are needed in order to continuously update the filter 
coefficient [3, 4, 5]. 

In order to investigate the ECG record of the patient in real-time, there is a chance that the ECG signal may be corrupted with 
noise. The main noise(artifact)available in the ECG contains: power-line interference (PLI), baseline wander (BW), muscle artifacts 
(MA) and motion artifacts (EM), that are generated by patient breathing, movement, power line interference, bad electrodes and 
improper electrode position. these artifacts strongly influences the ST section, degrades the sign first-rate, frequency resolution, 
produces huge amplitude indicators in ECG that can resemble PQRST waveforms and mask tiny capabilities which might be 
important for clinical tracking and diagnosis [1,2]. Most types of artifacts which affect ECG record are eliminated by band pass filters 
which may not give good result. The noise contaminated ECG records are acquired from the standard Massachusetts Institute of 
Technology-Beth Israel Hospital (MIT-BIH) data base. To enhance the accuracy and reliability for refine diagnosis the artefacts in ECG 
need to be decreased. Many methods have been implemented to eliminate the noise from noise contaminated ECG using static 
filters for example high pass, low pass and notch filters. However, the disadvantage of the fixed filter is that filtering or pre- 
processing of an ECG need to be done by understanding the noise cut off frequencies present in the signal. To overcome the 
restrictions of static filters, various methods have been developed. Adaptive filtering uses variety of algorithms such as LMS & NLMS 
algorithms etc. The performance execution of these algorithms is based on filter tap length (M) and the learning rate parameter (). 
The main drawback of using FIR filter in adaptive filtering is high computational load and slower convergence [6, 7]. 

Figure 1 defines the essential downside and noise cancelling answer. A signal s(n) is communicated over a channel to a detector 
that also collectsainterferencen1(n) unrelated with the signal. The first input to the ANC is noise contaminated signal d(n)i.e. 
s(n)+n1(n). The second detector accepts a noise signal n2(n) unrelated with the s(n)however correlative with n1(n) that supplies the 
reference input to the ANC that is filtered to produce the estimate signal y(n) that’s a detailed duplicate of reference input. The 
output signal y(n) is subtracted from signal d(n) to produce the error signal as shown in equation (1). This adaptive filter can be 
realized using various structures; the most frequently used structure is transversal finite impulse response (FIR). [8,9] 


e(n)=d(n)-y(n) (1) 
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Figure 1 Adaptive Noise Canceller 


Traditional filters such as adaptive filters [10], sign based normalised adaptive filters [11] were proposed in the literature to 
minimise artifacts. Different methodologies for ECG denoising include new variable step size NLMS [12,13] and EEMD-BLMS 
methods [14]. Promising performance are acquired by non linear adaptive algorithms [15,16], recently hybrid techniques have been 
proposed to eliminate noise from ECG signals using cascaded adaptive filters [17,18]. Filter banks (FB) decompose a digital signal 
into different frequency bands. On the basis of time frequency resolution, FB’s can be classified into two categories “e., uniform filter 
banks (UFB) and non uniform filter banks (NUFB). UFB provides fixed and uniform time frequency resolution. Whereas non uniform 
filer bank provides non-uniform and variable time frequency resolution lead to better performance and reduced arithmetic 
complexity in some applications like audio analysis, ECG signal enhancement [19,20]. Different methods of decomposing signals in 
to sub bands have become prominent and were proposed in the literature [21, 22]. 

A filter bank is used to decompose a discrete time signal into several sub bands by means of a set of analysis filters. Sub band 
adaptive filter (SAF) structures have been proposed to overcome these problems of adaptive filters. As a result, adaptive filtering 
using sub band coding becomes a good option for plenty of adaptive systems. Sub band adaptive filtering belongs to two fields of 
signal processing namely adaptive and multirate signal processing .Recently different noise cancellation methods are proposed 
using SAFs include variable step size sign SAF [23], Variable individual step size SAF [24], New normalised SAF [25]. An optimised 
cosine modulated NUFB design approach has been proposed by Kumar,A.,G K Singh et., all [26] .This proposed method uses FBs to 
decimate the input signal into a number of frequency bands, each serving as an unbiased input to SAF. Sub band signals are usually 
down sampled in a multirate system, which leads to a whitening of the input signals and therefore an improved convergence 
behaviour of the adaptive filter system is expected [5,39,40,41,42]. In conventional SAF’s each bound uses an individual adaptive sub 
filter in its own adaption that decreases the convergence rate. To solve these structured problems a multiband structured SAF 
(MSAF) are developed in which the full band adaptive filter's tap weight vectors are updated by a single adaptive algorithm using 
sub band signal [20,27]. The non-uniform filter bank SAF (NUFBSAF) is developed to achieve a better convergence performance by 
adapting the band width of analysis filter through proper selection of decimation factors. The objective of this paper is to develop 
non uniform multi band structured sub band adaptive filter which can improve the performance of the traditional ANC system, to 
analyse the application of SAF to the problem of noise cancellation problem in ECG. 


2. PROPOSED DESIGN METHODOLOGY 


LMS is most generally used adaptive algorithm than recursive least square algorithm and kalman filtering algorithm which makes 
use of a gradient vector to approximate a time-varying signal. A FB accommodates a set of analysis filters which decomposes the 
input signal of bandwidth into sub band signals. The sub bands give data from different frequency bands thus; it is conceivable to 
perform time and frequency dependent processing of the input signal. Since the sub bands are down sampled, processing is to be 
performed at a computationally efficient rate to analyze the ECG. A general objective is to create one arrangement of preprocessing 
filters which is useful in a noise cancellation tasks for ECG processing. In this manner, a sub band coding based adaptive algorithm 
calculation involves decomposing a signal into frequency sub bands and handling these sub bands as indicated by the current 
application. The concept of multiband structured-sub band adaptive filter (MSAF) is presented in this section.In the proposed design 
the primary input signal given to the first SAF consists of Noise contaminated ECG record denoted as d(n) in Figure.2. The secondary 
signal given to second SAF is noise signal n2(n) denoted as u(n). The full band input signal u(n), primary input signal or desired 
response signal d(n) and filter output signal y(n) are decomposed into N sub bands by means of analysis filters Hi(z); for (=0, 1,2...(N- 
7).In thisFigure.2 HoH7H2....Hy-7and FoF 7F2....Fy-7are analysis &synthesis filters of N channel perfect reconstruction (PR) filter bank 
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respectively. These sub band signals are decimated to a lower rate using same factor and are processed by individual sub band 
adaptive sub filters W;(z).[20] 
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Figure 2 Block diagram of a multiband structured sub band adaptive filter (adapted from Reference [20]) 


Non uniform filter banks have non uniform frequency partition. One method of building NUFB is to cascade UFB in a tree 
structure using a two channel FB as basic building blocks. The two kinds of tree-structured filter banks are regular and irregular. In a 
regular sub band tree-structure, sub bands in each level are decomposed by the same filter banks where as in irregular sub band 
tree-structure, sub bands in each level are decomposed by different filter banks or only some sub bands are decomposed. Irregular 
sub band tree structure is one possible solution for non uniform filter banks as mentioned above. The distinct advantage of the 
equivalent form of the tree-structured FB is the possibility of realizing analysis and synthesis filters for more important tasks than 
simple band separation [19]. As an example of this method, NUFB with decimation factors (16,16,8,4,2) is illustrated. The proposed 
tree-structured non uniform filter bank's analysis section is shown in Figure 3. 
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Figure 3 Proposed block diagram of a tree structured NUFB 
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The proposed structure makes use of low pass FIR prototype Parks-McClellan optimal equiripple filter design in both analysis 


and synthesis FB that give near-perfect reconstruction by permitting little measure of distortion at the output. The generalized 
structure of N-channel NUFB based on tree structured approach having decimation NoNz,N2---Ny-7 for each band then the 
decimation factors must fulfill the accompanying condition [21] 

N-1 1 _ 
what = 1 2) 
In N-channel NUFB the PR is possible if 
DidlH(e™)|° =1 foro<w <= (3) 
Where H;,(e/”) is frequency response of k * filter in equivalent NUFB parallel form 
For tree structured NUFB design having decimation factors (16,16,8,4,2) the PR condition can be achieved by using following 
equation 


|Ho(ei) |? + |H,(e)|? + |H2(e”)| + |H3(e™)|" 7 |Ha(e”)| =1for0<w == (4) 


Here Ho(z),H1(z), H2(z),H3(z),H4(z) analysis filters with the following relations [21] 


H(z) = Hy1 (2) H21 (27) 31 (2*) Hai (2°) (5) 
Hy (2) = Hy (2)H 21 (2?) 31 (24) Ha (2°) (6) 
H2(z) = H14(2)H21 (2”)H32(z*) (7) 
H3(z) = H44(z)H22(z*) (8) 
Hy4(Z) = Hy2() (9) 


Where H77, H21, H31, Haare low pass filters in the first second, third and fourth stages respectively and H12, H22, H32, H42 are high 
pass filters in the first second, third and fourth stages respectively. 


2.1. Analysis of conventional adaptive algorithms 

The traditional ANC shown in Figure.1incorporates filter whose input parameters are d(n)and u(n). Here u(n) is the time delayed 
input vector values, u(n) = [u(n),u(n-1),u(N-2),.........006 u(n-M +1)]7, the vector w(n)=[wo(n),w;(n),w2(n),.....wm-1(n)]™ represents 
the coefficients of the filter at time index n. The execution investigation of ordinary adaptive filter algorithms requires five particular 
computational strides in every emphasis as takes after.[5] 


1. The estimated response y(n)of the FIR filter calculated as 
y(n) = Yin=ow (nju(n — m) = w"(n)u(n) (10) 


2. The difference between primary input signal and estimated response signal of the adaptive filter £e. error signal is 
e(n) = d(n) — y(n) (11) 


3. The minimization of the mean-square error (MSE) ie. Performance function, which is defined as 


J = Efe?(n)} (12) 


4. For a given filter coefficient w(n) with reference input signal u(n) and primary input signal d(n) , the MSE can be calculated as 
J = E{d?(n)} — 2p’w + w" Rw (13) 


Where R is the autocorrelation function of input signal calculated as 
R= E {u(n) u"(n)} (14) 
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And p is the cross-correlation array between the input signal vector and the primary input signal calculated as 
p= E{d(n)u(n)} (15) 


5. The filter coefficients or weight vectors of the filter are changed for the next iteration for LMS adaptive algorithm by using 
equations 
w(n +1) = w(n) + pe(n)u(n) (16) 


Here yu is chosen such that to guarantee the stability of the algorithm in the range 
(17) 


0<u< oP, 


Here Py is average power of the reference signal u(n) calculated as 
P, =u" (n)u(n) (18) 


2.2. Analysis of MSAF algorithms 

Figure.2 shows two equivalent structures for the selected input portion of the multiband system as indicated that the full band 
adaptive filter W(z)of length ™ is stationary (i.e. the adaptation is frozen), it can be transposed with the analysis filter bank as each 
sub band signal occupies only a portion of the original frequency range. The N sub band signals of MSAF algorithms 
requireparticular computational strides in every emphasis as takes after.[20] 


It is important to note that n refers to the time index of original sequence and k denotes the time index of decimated signal. 


1. The adaptive filtered output signal is 

yin) = Linz w (n)u(n —m) =w"(n)uj(n) (19) 

FOF THO Ay os cetices N-1 

2. The adaptive filtered decimated output signal is 


Viv lk) = yi(KN) = Linzo w (k)u(k — m)=w" (k)u;(k) 

for tS 0A evens: N-1 (20) 

Where w(k) is the weight vector of the full band adaptive filter and ui(k) is the regression (signal) vector for the it” sub band at every 
time instant k, each sub band signal vector u;(k) is packed with N new samples and (M —N) old samples to produce an output 


sample yi p(k). 


3. The shift register performs both decimation and packing operations to form the sub band vectors defined in equation (21) 
uj(k) = [uj(KN), uj(KN — 1), uj(KN — 2)......u,(KN —M +1)]? (21) 


4. The sub band decimated error signal e;p(k) is the difference between the filter output y; p(k) and the desired response 
dj p(k)expressed as 

e,0(k)=di,p(k) — Yip (k)=dip(k) — w" (k)uj(k) (22) 

FOF T= Oy oscciice: N-1 

5: The error signal vector for N sub bands, 

ep(k) = [eon(k), e1,n(k), €2,p C3 Poernerers ey—1,p(k)]" (23) 

6. The error signal can be expressed in a compact form as 

ep(k) = dp(k) — U"(k)w(k) (24) 


wheredp(k) is N X 7 decimated desired response vector and U(k) is M X N sub band signal matrix defined as follows 
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U(k) = [uo (k), uy (K), Ug (K),.. cece uy_—1(k)] (25) 


dp(k) = [do (k), di,p(k), donk), .--e- eee dy-1,p(k)]" (26) 


The multi band MSE functions as a performance measure for the MSAF algorithm. By means of the Wiener filter theory, the 
multiband MSE function can be shown to be equivalent to the classical full band MSE function. Hence, use of the multiband MSE 
function allows a fair comparison of the MSAF algorithm with other adaptive algorithms that are measured by the MSE function. 


i The multiband MSE function is defined as the average of the mean-squared values of the sub band estimation errors 
expressed as [20] 

In= rit Efe; p7(n)} = og? — 2p’w + w? Rw (27) 

8. The multiband MSE function can be written in the following parameters variance (o42) of the primary signal auto 


correlation matrix (R) and cross correlation matrix (p) expanded form: 


oq? = E{d?(n)} = —F{lldp(k)II7} (28) 
p = E{u(n)d(n)} = -E{U(k)dp(&)} (29) 
R= Efu(n)u" (n)} = SE{U(KUT (k)} (30) 


A recursive weight-control mechanism for the MSAF from one iteration to the next iteration, coefficients of filter should be 
changed in a minimal manner. This basic principle is used to design a constraint optimization criterion for deriving the adaptive 
algorithm. In the MSAF shown in Figure.2 with the sub band signals sets {U(k),dp(k)},a criterion that ensures convergence to the 
optimum solution after a adequate number of iterations is to have the updated tap weight vector w(k+ 17) using LMS Algorithm in all 
N sub bands at each iteration k as follows 


d., 
wk + 1) =w(k) +h Td ieaa Ci (*) (31) 
Where uy is learning rate parameter and it should be selected in the stability bound ie equation (15) 
Where a is a small positive constant used to avoid possible division by zero. 


Proposed NUFB-MSAF algorithm summary 
FON 'S= 1) Qoecccesdcans N where s is stages 


Analysis filters Hi(z) 
for (=0,1,2........ N-7 


Hig=Hale) 2). ce) 
For k=O, 1, 2y.cccscee kN where kN=n 


Error estimation: 
€p(k) = dp(k) — UT (k)w(k) 


Normalization matrix: 
A(k) = diag[U" (k)U(k) + a] 
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Tap-weight adaptation: 
wk +1) = wk + 1) + WU (K)A 1 (hep (k) 


Input signal Band partitioning: 
UT (k) = HT™A(KN) 


Desired signal Band partitioning: 
dp(k) = H™A(KN) 


Synthesis: 
e(kN) =Fep(k) 


Parameters: 

| &h - Basic Low pass filter& high pass filter notation 
M - Number of adaptive tap weights 

N - Number of sub bands 

L - Length of the analysis and synthesis filters 


Variables: 

UT (k) = [Uf (k), UZ (k — 1)). 

UZ(k —1) - first M-N columns of U7 (k — 1) 

A(KN) = [a(KN), a(KN — 1),a... ....a(CKN — N — 1)] 

a(kN) = [u(KN), u(kKN — 1),u(KN — 2),..u(kKN —L —1)]". 
d(kN) = [d(KN),d(kKN — 1),d(kN — 2)........d(kKN —L—1)]’. 


3. RESULTS AND DISCUSSION 
In this simulation the benchmark MIT-BIH arrhythmia database [28] was utilized to test the execution of various UFB & NUFB 
structured MSAF's adaptive algorithms for ECG denoising and it comprise of 48 half hour excerpts of two channel ambulatory ECG 
recordings, which were acquired from 47 subjects, including 25 men aged 32-89 years, and women aged 23-89 years. The 
recordings were sampled at 360 samples per second per channel with 11-bit resolution over a 10 mV range. The simulations were 
done by collecting 3600 samples of ECG recordings. In this simulation the proposed methodology has been implemented using 
Parks-McClellan optimal equiripple linear phase (real, symmetric coefficients) low pass FIR filter design with length 32,cut off 
frequency 100 Hz and sampling frequency 360 Hz. A data set of five ECG records: data100, data105, data108, data203 and data228 
are considered to guarantee the consistency of the results. The secondary noise signal n2(n) shown in Figure.1 is taken from noise 
generator. A synthetic PLI with Imv amplitude is simulated for PLI cancellation. In order to test the filtering capability in non- 
stationary environment real PLI noises with 3600 samples are considered. These are taken from MIT-BIH normal sinus rhythm 
database (NSTDB) [29]. 

The proposed UFB and NUFB structured MSAF-ANC are compared quantitatively by quality assessment parameters Mean Square 
Error (MSE), Root Mean Square Error (RMSE), Signal-to-Noise Ratio Before Filtering (SNRBF), Signal-to-Noise Ratio After Filtering 
(SNRAF)and distortion are outlined in Table.1. 


Table 1 Proposed design's quality assessment parameters 


PARAMETERS FORMULA 
N-1 

y. [s(n) — e(n)]? 
7 


0 


1 
Mean Square Error 7 


Root M.S.E 


N- 
i> [s(n) — e(n)]? 
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N-ALs(n)P 

10 logio Sttaq@) — SCO 
N=ALs(n)? 

10 logio 1 


N-1N-1 


10 logy y y [s(n) — e(n)]? 


n=0 n=0 


SNRBF ) 


SNRAF 


Distortion 


3.1. Power-line Interference (PLI) Cancellation using MSAF-LMS algorithm 

In this proposed design 3600 samples of clean ECG signal records from MIT-BIH data base which are corrupted with power line 
interference noise with frequency60 Hz and the sampling frequency of 360Hz is applied as primary input signal to the ANC shown in 
Figure 1. The following figures 4-15show the PLI Noise elimination. The NUFB structured MSAF’s achieve good SNR; MSE and RMSE 
values over UFB structured MSAF's adaptive algorithms. Table’s 3-6 presents ‘quantitative parameter values for five record numbers 
from MIT-BIH data base. The simulation results state that NUFB structured MSAF's adaptive algorithms has greater efficiency than 
UFB structured MSAF’s adaptive algorithms. These results for MIT-BIH data base record number 105 are shown in Figures 4-15. 


Table 2 Simulated SNR values in dB 


SNR Values 100 105 108 203 228 Avg.Val 

Before Filtering 0.2144 0.712 -0.0059 2.9671 -1.0532 0.5669 

Three Channel(3,3,3)UFB 11.5386 11.8879 11.2808 13.7657 10.5938 11.81336 
Four Channel(4,4,4,4)UFB 11.9735 12.4391 11.7256 14.0723 11.2329 12.28868 
Five Channel(5,5,5,5,5)UFB 11.1846 11.6595 10.9155 13.5753 10.2149 11.50996 
Three Channel(4,4,2)NUFB 14.6593 15.0662 14.4046 16.71 13.9338 14.95478 
Four Channel(8,8,4,2) NUFB 15.7315 16.9193 15.4913 18.9523 14.8402 16.38692 
Five Channel(16,16,8,4,2)NUFB 21.791 21.214 21.833 22.9534 20.7571 21.7097 


The performance of UFB structured MSAF’s and NUFB structured MSAF’susing LMS algorithms in terms of SNR are shown in 


Table.2. The simulation results state that NUFB structured MSAF algorithms has greater efficiency than UFB structured MSAF 
algorithms. As shown in Table.2 three channels UFB has on average SNR of 12.44774dB for 5 records and three channels NUFB has 
on average SNR of 14.95478dB. Similarly for four channel UFB has on average SNR of 12.28868dB for 5 records and four channels 
NUFB has on average SNR is 16.38692dB, for five channel UFB has on average SNR is 11.50996dB for 5 records and five channel 
NUFB has on average SNR is 21.7097dB.It is clear from Table.2 that NUFB structured MSAF’s adaptive algorithms outperform UFB 
structured MSAF algorithms in approximating the ECG noises. 


Table 3 Simulated MSE values 


Mean Square Error Values 100 105 108 203 228 Avg.Val 
Three Channel(3,3,3)UFB 0.0584 0.1831 0.0557 0.457 0.0672 0.16428 
Four Channel(4,4,4,4)UFB 0.062 0.186 0.0579 0.4493 0.0685 0.16474 
Five Channel(5,5,5,5,5)UFB 0.0574 0.1825 0.0546 0.4515 0.0653 0.16226 
Three Channel(4,4,2)NUFB 0.0884 0.2569 0.082 0.6152 0.0962 0.22774 
Four Channel(8,8,4,2) NUFB 0.1631 0.4716 0.1573 1.038 0.1694 0.39988 
Five Channel(16,16,8,4,2)NUFB 1.0749 1.2248 1.1383 2.414 0.9461 1.35962 


Table 3 shows the performance of UFB structured MSAF and NUFB structured MSAF using LMS algorithms in terms of MSE. It is 
clear from Table.3 that UFB structured MSAF adaptive algorithms have minimum MSE values than NUFB structured MSAF adaptive 


algorithms. 
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Table 4 Simulated RMSE values 


Root Mean Square ErrorValues 100 105 108 203 228 Avg.Val 
Three Channel(3,3,3)UFB 0.2416 0.4279 0.2361 0.676 0.2591 0.36814 
Four Channel(4,4,4,4)UFB 0.249 0.4313 0.2406 0.6703 0.2616 0.37056 
Five Channel(5,5,5,5,5)UFB 0.2396 0.4272 0.2336 0.6719 0.2555 0.36556 
Three Channel(4,4,2)NUFB 0.2973 0.5068 0.2863 0.7843 0.3101 0.43696 
Four Channel(8,8,4,2) NUFB 0.4039 0.6867 0.3966 1.0188 0.4116 0.58352 
Five Channel(16,16,8,4,2)NUFB 1.0368 1.1067 1.0669 1.5537 0.9727 1.14736 


Table 4 shows the simulated MSE values of UFB structured MSAF and NUFB structured MSAF using LMS algorithms. It is clear 
from Table 4 that UFB structured MSAF algorithms have minimum RMSE values than NUFB structured MSAF algorithms. 


Table 5 Simulated distortion values 


Distortion Values 100 105 108 203 228 Avg.Val 
Three Channel(3,3,3)UFB -12.3381 -7.3742 -12.5378 -3.4006 -11.7291 -9.47596 
Four Channel(4,4,4,4)UFB -12.0774 -7.3041 -12.3732 -3.4744 -11.646 -9.37502 
Five Channel(5,5,5,5,5)UFB -12.4117 -7.3867 -12.6298 -3.4533 -11.8519 -9.54668 
Three Channel(4,4,2) NUFB -10.5349 -5.9027 -10.8634 -2.1099 -10.1704 -7.91626 
Four Channel(8,8,4,2)NUFB -7.8743 -3.2645 -8.0329 0.1621 -7.7096 -5.34384 
Five Channel(16,16,8,4,2) NUFB 0.3137 0.8808 0.5625 3.8274 -0.2404 1.0688 


Table 5 shows the performance of UFB structured MSAF and NUFB structured MSAF using LMS algorithms in terms of distortion. 
It is clear from Table 5 that UFB structured MSAF adaptive algorithms have minimum distortion values than NUFB structured MSAF 
algorithms. 
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Figure 4 Simulation results using MSAF-LMS algorithm for record number 105 (i) ECG record from MIT-BIH database (red) (ii) PLI 
contaminated ECG record (magenta) (iii) Sub band adaptive filtered signal using three channel (3,3,3) uniform filter bank MSAF 
using LMS algorithm (red) 
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Figure 5 Three channel UFB with decimation factors (3,3,3) simulated results in dB (i) Magnitude spectrum of analysis FB (ii) 
Magnitude spectrum of synthesis FB (iii) Amplitude variation of PRE (iv) Reconstruction error Amplitude 


The simulation results of a proposed method with proto type filter for three channel UFB with decimation factors (3,3,3) is shown 
in Fig.5.The peak reconstruction error (PRE) is reduced appreciably. The average PRE obtained in three channel UFB is 0.08169 dB. 
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Figure 6 Simulation results using MSAF-LMS algorithm for record number 105 (i) ECG record from MIT-BIH database (red) (ii) PLI 
contaminated ECG record (magenta) (iii) Sub band adaptive filtered signal using four channel (4,4,4,4) uniform filter bank MSAF 
using LMS algorithm (red) 
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Figure 7 Four channel UFB with decimation factors (4,4,4,4) simulated results in in dB (i) Magnitude spectrum of analysis FB (ii) 


Magnitude spectrum of synthesis FB (iii) Amplitude variation of PRE (iv) Reconstruction error Amplitude 


The simulation results of a proposed method with proto type filter for four channel UFB with decimation factors (4,4,4,4) is 
shown in Fig.7. The peak reconstruction 
0.0414cB. 


error (PRE) is reduced appreciably. The average PRE obtained in four channel UFB is 
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Figure 8 Simulation results using MSAF-LMS algorithm for record number 105 (i) ECG record from MIT-BIH database (red) (ii) PLI 


contaminated ECG record (magenta) (iii) Sub band adaptive filtered signal using five channel (5,5,5,5,5) uniform filter bank MSAF 
using LMS algorithm (red) 
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Figure 9 Five channel UFB with decimation factors (5,5,5,5,5) simulated results in dB (i) Magnitude spectrum of analysis FB (ii) 
Magnitude spectrum of synthesis FB (iii) Amplitude variation of PRE (iv) Reconstruction error Amplitude 


The simulation results of a proposed method with proto type filter for five channel UFB with decimation factors (5,5,5,5,5) is 
shown in Fig.9. The peak reconstruction error (PRE) is reduced appreciably. The average PRE obtained in five channel UFB is 
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Figure 10 Simulation results using MSAF-LMS algorithm for record number 105 (i) ECG record from MIT-BIH database (red) (ii) PLI 
contaminated ECG record (magenta) (iii) Sub band adaptive filtered signal using three channel (4,4,2) non uniform filter bank MSAF 


using LMS algorithm (red) 


pase | 9 Z 


ARTICLE 


(a) Amplitude responses of the analysis filters of filter bank in dB 
T T T T T 7 T 


Magnitude in dB 


Magnitude in dB 


20 


of filter bank in dB 


T T 20 


Magnitude in dB 


(b) Amplitude responses of the synthesis filters 
T 7 T T T T 


T T 


(d) Variation of reco: 
T T 


0.4 0.5 0.6 0.7 


Normalized Frequency 


nstruction error in dB 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0 0.1 
Normalized Frequency 
a (c) Variation of peak reconstruction error in dB 5 
T T T T T T T 
-20 1 
o 
3 
£e 
o 
-40 5 4 3 
2 
& 
= 
-60 F | 
-80 f L L f f i f f ! -70 
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0 0.1 


Normalized Frequency 


fl 
0.4 0.5 0.6 0.7 
Normalized Frequency 


Figure 11 Three channel NUFB with decimation factors (4,4,2) simulated results in dB (i) Magnitude spectrum of analysis FB (ii) 
Magnitude spectrum of synthesis FB (iii) Amplitude variation of PRE (iv) Reconstruction error Amplitude 


The simulation results of a proposed method with proto type filter for three channel NUFB with decimation factors (4,4,2) is 
shown in Fig.11.The peak reconstruction error (PRE) is reduced appreciably. The average PRE obtained in three channel NUFB is 


0.03583cB. 
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Figure 12 Simulation results using MSAF-LMS algorithm for record number 105 (i) ECG record from MIT-BIH database (red) (ii) PLI 
contaminated ECG record (magenta) (iii) Sub band adaptive filtered signal using four channel (8,8,4,2) non uniform filter bank 
MSAF using LMS algorithm (red) 
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Figure 13 Four channel NUFB with decimation factors (8,8,4,2) simulated results in dB (i) Magnitude spectrum of analysis FB (ii) 
Magnitude spectrum of synthesis FB (iii) Amplitude variation of PRE (iv) Reconstruction error Amplitude 


The simulation results of a proposed method with proto type filter for four channel NUFB with decimation factors (8,8,4,2) is 


shown in Fig.13.The peak reconstruction error (PRE) is reduced appreciably. The average PRE 


0.023393cB. 
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Figure 14 Simulation results using MSAF-LMS algorithm for record number 105 (i) ECG record from MIT-BIH database (red) (ii) PLI 
contaminated ECG record (magenta) (iii) Sub band adaptive filtered signal using five channel (16,16,8,4,2) non uniform filter bank 
MSAF using LMS algorithm (red) 
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Figure 15 Five channel NUFB with decimation factors (16,16,8,4,2) simulated results in dB (i) Magnitude spectrum of analysis FB (ii) 
Magnitude spectrum of synthesis FB (iii) Amplitude variation of PRE (iv) Reconstruction error Amplitude 


The simulation results of a proposed method with proto type filter for five channel NUFB with decimation factors (16,16,8,4,2) is 
shown in Fig.11. The peak reconstruction error (PRE) is reduced appreciably. The average PRE obtained in five channel NUFB is 
0.0288cB. 

The simulation results of a proposed method using LMS adaptive algorithm with proto type filter for three channel, four channel 
and five channel UFB’s and NUFB’s are shown in above Figures. The peak reconstruction error (PRE) is reduced appreciably. The 
average PRE obtained for three channel, four channel and five channel UFB’s are 0.08169 dB, 0.0414 dB and 0.0454 dB respectively. 
The average PRE obtained for three channel, four channel and five channel NUFB's are 0.03583 dB, 0.0238393 dB and 0.0288 dB 
respectively. The performance comparison with published works is shown in Table 10. 


Table 10 Comparison of existing works 


SNO Author Year Method SNR in dB 
1 Gowri, Thumbur et.,al [14] 2014 | E- VSSLMS algorithm 12.3451 
Keergaard, Kevin et.,al [13] 2016 | EEMD-BLMS algorithm 13.1900 
Proposed design 2018 | NUFB-MSAF LMS algorithm 21.7097 


4. CONCLUSIONS 


This research paper presents the new comprehensive ANC system of ECG signals with robustness based on UFB &NUFB structured 
MSAF's using LMS algorithm. The proposed model is potentially a new realization structure form of ANC which guaranteed a more 
stable transformation in response to variants in input signal power. The theoretical analysis of NUFB structured MSAF system is 
carried out and simulations are performed using MATLAB. In order to analyze the performance of the proposed design a 
comparison has been made between six different ECG denoising schemes ie three channel, four channel and five channel UFB & 
NUFB structured MSAFs using LMS algorithm respectively. The proposed NUFB structured MSAF system consists of reference signal 
and primary signal as input parameters for which adaptive filtered estimated signal, error signal and coefficients of the filter are 
obtained as output parameters. The performance of proposed system is compared quantitatively by parameters SNR, MSE, RMSE 
and Distortion. Better filtering performance results are obtained by NUFB structured MSAF using LMS algorithm and also this 
algorithm guarantee the better estimation of noise. Computer simulation demonstrated that the proposed system gives improved 
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performance and achieves good adaptation. The SNR for various NUFB structured MSAF's was found to be higher than the UFB 
structured MSAF's. The five channel NUFB structured MSAF performs better SNR values than UFB structured MSAF using LMS 
algorithm. 
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